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EDITORIAL 


U 


« 


Chers co-propriétaires de PPC, 


Retenez dès à présent la date du 7 janvier 1989. C’est en effet la date de notre 
prochaine Assemblée Générale. Elle aura lieu comme d'habitude au Centre Jean 
Verdier, de 14 à 16h en salle 303. 


A ce propos, je tiens à vous rappeler qu’une seule personne s’est présentée pour 
faire partie du nouveau Bureau. C’est insuffisant. Avez-vous pensé à toutes les 
implications ? Pour ma part, j'en vois une qui ne me réjouis guère... Lisez-donc 
Particle en page 2. 


La Bibliothèque de Genève progresse considérablement, grâce aux efforts de 
Thomas Rouyer. Tous les possesseurs de HP-41 sont vivement invités à lire 
l’article dans les pages générales. 


Nous avons vu récemment le HP-42S, une nouvelle calculatrice compatible 
HP-41. Le format est celui d’une HP-27S. Le HP-42S reprend visiblement le coté 
calculatrice du HP-41 sans en reprendre le coté ordinateur. Une sortie 
imprimante infra-rouge, mais toujours pas de dispositif de sauvegarde de 
programmes ou autres entrées/sorties. Dommage. 


Enfin, je tiens à vous signaler que ce JPC fait 32 pages. C’est peu ? Oui. Il ne 


tient qu’à vous d’y remédier. Il nous reste quelques articles pour le mois 
prochain, mais ce n’est pas suffisant. Compris ? Alors, à vos plumes... 


Mais je parle, je parle, et personne ne m'’arrête. Je me tais donc, et vous laisse 
découvrir votre Journal. 


En espérant vous apporter une très bonne nouvelle le mois prochain... 


Pierre David (37) 





vous” 


qu'ils clisatent 


JPC 59 Page 1 


L'AVENIR DU CLUB 


Nous sommes le 21 octobre et, malgré les appels 
lancés à maintes reprises, une seule personne s’est 
présentée pour nous remplacer et faire partie du 
nouveau Bureau. 


Cette situation est très préoccupante. 


Avez-vous pensé à ce que signifie votre absence ? Si 
personne ne se présente, le Club sera obligé de cesser 
ses activités. Plus de Club, plus de JPC, plus de 
programmes, plus d'informations, plus de réunions. 
La fin. 


Or cela serait grave. Beaucoup de travail a été 
accompli. Que deviendra-t-il ? Le Club répond à un 
besoin, je m’en rends compte tous les joufs. 
Comment ferez-vous si le Club disparaît ? 


Comment se fait-il que personne ne se manifeste ? 
Toutes les personnes interrogées disent ne pas avoir 
le temps. Croyez-vous vraiment que nous ayons plus 
de temps libre que vous ? 


Enfin, nous nous posons une question de fond : 
tenez-vous au Club ? 


Si la réponse est négative, si personne ne se présente 
avant le 15 novembre dernier délai, il faudra bien se 


rendre à l'évidence et engager le processus de 
cessation des activités. 


Nous avons passé beaucoup de temps à nous occuper 
du Club. Le jour où il faudra aller à la Préfecture 
pour fermer le Club sera très pénible. Mais nous 
irons la conscience tranquille. Nous aurons fait tout 
notre possible pour ne pas en arriver là. Mais s’il le 
faut... 


Réagissez ! Il est encore temps. 


Olivier Arbey (118) 
Jacques Baudier (192) 
Pierre David (37) 
Jean-Jacques Dhénin (177) 
Eric Gengoux (108) 
Janick Taillandier (246) 
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LA BIBLIOTHEQUE DE GENEVE 


Vous le savez depuis le mois précédent, la 
Bibliothèque de Genève a été reprise par Thomas 
Rouyer. Qu'est-ce que la Bibliothèque de Genève ? 
Eh bien, il s’agit d’un vaste ensemble de programmes 
pour HP-41, datant des temps Zzéroïques, et 
maintenant en notre possession. 


Mais il faut s’en occuper. Et Thomas a déjà trié les 
programmes en plusieurs séries : 

- analyse numérique - mathématiques 
- géométrie et trigonométrie 

- génie civil 

- antennes 

- astronomie - optique 

- électricité - électrotechnique 

- médecine 

- industrie pétrolière 

- théorie des nombres 

- jeux 

- programmes d'informatique 

- programmes financiers - comptabilité 
- navigation aérienne et marine 

- mécanique - physique 

- etc. 


Si vous êtes intéressés, faites nous le savoir. Il serait 
dommage de ne pas profiter de ce travail. 


Par ailleurs, Thomas cherche des pochettes plastiques 
(format A4) avec logement pour cartes magnétiques 
afin de classer ces programmes et leur 
documentation. Si vous savez comment vous en 
procurer, Thomas pourra continuer son travail. 


Avis donc à tous les possesseurs de HP-41. Vous avez 
là une gigantesque réserve avec un travail mâché, 
profitez-en.. 


Pierre David (37) 


Pierre Picheret 
Calmon Aiguefonde 
81200 Mazamet 

Tel : 63 61 27 11 


Vend : 


Imprimante HP-IL HP82905B : 1000 FF, imprimante 
ThinkJet HP2225B : 1500 FF. 


Janick Taillandier 
335 rue Lecourbe 
75015 Paris 

Tel : (1) 45 54 18 83 


Vend : 
Module  HHP 
HP-71 : 1000 FF. 


32Ko Ram / 32Ko Eprom pour 


François Laissus 

14 avenue Jean Moulin 
75014 Paris 

Tel : (1) 45 39 06 31 


Vend : 

Lecteur de cassettes HP82161A : 1000 FF, module 
HP-IL pour HP-41 en prime. Convertisseur HP-IL 
HP82166 (jamais servi) : 500 FF. 





BILLET DOUX 
(ACTE Ill) 


Ils sont pas possibles, à la SNCF ! Ils ont encore 
augmenté leurs tarifs grandes lignes le premier juin. 
Alors, je suis obligé de mettre à jour mes 
programmes pour la HP-41 et le HP-71 (voir JPC 36 


et JPC 44)... 


Profitons-en pour vous en proposer une version pour 
la HP-27S avec, grâce à sa fonction SOLVE, deux 
nouvelles possibilités : 


- répondre à la question suivante : quelle réduction 
faudrait-il, toutes choses égales par ailleurs, pour 
sortir à tel niveau de prix (par exemple, le coût du 
même trajet, effectué en voiture ou en autocar, ou 
même en avion...) 


- avoir une idée, connaissant le prix du billet, sa classe 
et sa réduction, de la distance ayant servi à l’établir.. 


Les lecteurs de mes deux précédents articles se 
souviennent probablement de la complexité de la 
formule de calcul du prix, qui comportait quatre 
grandes particularités : 


- détermination d’une distance dite palière : médiane 
de la tranche de 5 ou 10 km dans laquelle tombe la 
distance réelle, 


- choix, selon distance, de l’une des trois formules 
possibles la dernière étant, au-delà de 350 km, 
nettement dégressive.., 


-limitation, aux courtes distances, du poids de 
l'élément fixe dans le prix final du billet, 

-enfin, formules d’arrondi et de minimum de 
perception, assez spéciales. 


Pour ne pas trop compliquer les choses, on limitera la 
précision du calcul au Franc près, ce qui suffit pour 
un calcul courant, disons pour un trajet de plus de 
30 km. et on fera l'impasse sur certains détails : 
distance palière (il vous faudra la calculer 
approximativement(")), et arrondi précis (tout sera 
fait au Franc immédiatement supérieur). 


Dès lors, on va pouvoir utiliser à plein les possibilités 
du menu f[SOLVE], et notamment les clauses IF ; 
l'expression résultante n’est pas trop lourde, 
jugez-en : 
BIL = 0.01x(100-RED)x(2.5*CL)x... 
...IF(DIST<300:DISTx0.4572+10.30:... 
+... IF(DIST>349:DISTx0.37+39.00:... 
-..DISTx0.4212+21.00)) 


ce qui donne un menu très simple : 
(BIL] (RED] [CL] [DIST] 


On peut apporter une petite amélioration pour sortir 
le résultat en Francs entiers (mais ce n’est pas 
fondamental), en écrivant : 


BIL = 1P(0.99+0.01x(100-RED)x(2.5*CL)x... 
-..IF(DIST<300:DISTx0.4572+10.30:... 
...IF(DIST>349:DISTx0.37+39.00:... 
-..DISTx0.4212+21.00))) 


ce qui marche quel que soit le mode d’affichage 
utilisé, à la différence de la fonction RND du menu 
PARTS... 


Quelques explications sur les coefficients utilisés, qui 
vont permettre de mettre à jour au passage les 
programmes déjà publiés pour HP-41 et HP-71 (voir 
supra) : 
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| Taxe dite | Taxe dite | 

| complémentaire | kilométrique | 
mms tonne mm mm mme 
de 1 à 299 km | 10.30 F | 0.4572F | 
de 300 à 349 km | 21.00 F | 0.4212F | 
350 kmet plus | 39.00 F | 0.3700F | 
Sésame PER mer dpt SES ES dec LUS 
Les éléments ci-dessus concernent la deuxième classe, 
les prix en première étant toujours égaux à 1,5 fois 
ceux de seconde. 


Pour la HP-41, seul le programme PBASE devra être 
modifié (il charge les prix de base et seuils de 
déclenchement dans les registres RO1 à R12 de la 
HP-41). Aucun des 3 autres programmes, SNCF, D1 ou 
T2, n’est donc à modifier (voir JPC 44 page 8). 


Pour le HP-71, seules les lignes 90 et 100 du 
programme unique (j'allais dire commun...) sont 
affectées : 


90 PO(1)=10.30 à PO(2)=21 à PO(3)=39 
100 P1(1)=0.4572 à P1(2)=0.4212 à P1(3)=0.37 


Pour la bonne règle, modifier le nom du 
sous-programme, qui  s’appellera désormais 
SUB TAX688(C,R,B,D,P), et la remarque en ligne 30, le 
niveau tarifaire étant celui du premier juin 1988. 


Voilà. Un dernier mot quand-même : je ne pensais 
pas, quand j'ai écrit ce petit programme, que le menu 
SOLVE de la HP-27S était aussi souple et rapide (pas 
plus d’une seconde...). Un vrai plaisir ! 


Eric Gengoux (108) 





(*) Rappel : jusqu’à 199 km, les paliers sont de S km, la distance 
médiane se termine par un 2 ou un 7 ; au-delà, les paliers sont de 
10 km, leur distance supérieure se termine par un 9, donc la 
distance médiane par un 4. Exemple : Paris-Lyon, 519 km, utiliser 
514 km. — 


APE EEE STE Te. 
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HP28 


P. Heïlbronn Choc en retour 6 
D.T. Khoi-Nguyen Définitions de caractères 6 
S. Vaudenay Assembleur sur HP-28C 1CC 7 
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CHOC EN RETOUR 


Trois coquilles se sont glissées dans JPC 57, dans 
l’article Ajustement de courbes, en page 7 : 

- première colonne en sixième ligne, 

au lieu de {mR}, lire ©). 


- première colonne en deuxième ligne avant le bas : il 
manque une apostrophe (") à droite de la flèche + 
(ou, si l’on veut l'expliquer autrement, à gauche de 
IFTE). 


- deuxième colonne en onzième ligne : le dernier mot 
de la ligne n’est pas STO mais STOS. 


Dans JPC 57, page 10, rubrique Chocs en retour. Je 
remercie Paul Courbis de ses utiles indications basées 
sur les remarques dans l'article Le facteur est constant 
paru dans JPC 56. Je suis en plein accord avec la 
première indication. En revanche, j'ai essayé le 
programme ASSEMBLEUR sur une HP-28S avec l’adresse 
indiquée (#20000) et je n’ai pas réussi à implanter les 
libellés des fonctions arithmétiques comme sur une 
HP-28C. J'ai subi un Memory Lost. Ai-je commis une 
erreur de programme ou de manipulation ou bien 
plusieurs versions de HP-28S sont-elles sur le 
marché ? 


Philippe Heïlbronn (233) 





DEFINITION DE CARACTERES 


Je suis lecteur de JPC où je trouve beaucoup de trucs 
pour mon calculateur HP-28C. Bien qu’étant en 
dehors du Club, j'aimerais contribuer à l’élaboration 
du Journal en vous adressant ce programme de 
création de caractères selon la méthode employée par 
le HP-85B avec la fonction BPLOT$. 


« D! STO 
‘F! STO 
D SIZE + s 
« CLLCD 
15 
FOR Y 
D Y GET R-+B -+STR 
DUP SIZE 3 SWAP SUB DUP SIZE + ct 
«it 
FOR X 
c X X SUB X 
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IF 1 == 
THEN DROP 
ELSE 
IF mn == 
THEN 
X.I*FRE+SY-.1* 
F IM + R+C PIXEL 
END 
END 
NEXT 
» 
NEXT 

» 

{ D F PPAR } PURGE 
» 
ICHR' STO (mode tBIN]) 
Le principe du programme est simple : le caractère 
est divisé en lignes et chaque ligne remplie de chiffre 
1 et/ou 0 (1 pour une casse pleine et 0 pour une vide) 
constitue un nombre binaire. 


Utilisation : 

1) Reporter le dessin du caractère sur un quadrillage, 
2) Remplir les cases avec 1 si la case est pleine et 0 si 
elle est vide, 

3) A gauche du dessin, accoler une colonne de 1 qui 
sert de repère lors du tracé (mais n’apparaît pas), 

4) Transformer chaque nombre binaire en un réel par 
t8-R] et les placer dans une liste, 

5) Entrer un nombre complexe pour centrer le dessin. 
Le pixel correspondant se trouvera au coin gauche en 
bas, et 

6) Faire [SWAP] puis [CHR]. 


Exemple à cCHR] : un sigle bien connu... 


1001111101110 5102 
1011111011100 6108 
1111111011100 8156 
1111110011100 8092 
1111100111000 7992 
1111000111000 7736 
1111000111000 7736 
1111001111111 7807 
1110001111111 7295 
1110001110111 7287 
1110011101110 7406 
1110011101110 7406 
1110011101110 7406 
1110111011101 7645 
1111000000001 7681 
1111000000001 7681 
1111100000011 7939 
1111110000011 8067 
1111111000011 8131 


1011111110111 6135 
1001111110111 5111 pi 


1111111111100 8188 
1001110000000 4992 
1000001111111 4223 
1000000100000 4128 
1000000011111 4127 
1000000001000 4104 
1000000001111 4111 
1011111110100 6132 
1011111110111 6135 
1011101110100 6004 
1111011100111 7911 
1111011100100 7908 
1111111100111 8167 
1111111000100 8132 
1110000001111 7183 
1110000001000 7176 
1100000011111 6175 
1100000100000 6176 
1100001111111 6271 
1001110000000 4992 
1011111111100 6140 p2 


« CLLCD (-0,5;-1) P1 CHR (0,7;-1) P2 CHR » [EVAL] 


Note : éliminer CLLCD dans (CHR. 


Dang-Trung Khoi-Nguyen 





ASSEMBLEUR SUR HP-28C 1CC 


Aux lésés.. 


Vous devez être surpris de me retrouver dans ces 
colonnes, mais rassurez-vous, je reviendrai sévir dans 
les colonnes réservées au HP-71 dans peu de temps. 
Si je suis venu ici, c’est pour redonner espoir aux lésés 
de la version 1CC, jusqu'alors sous la tyrannie des 
1BB et des HP-28$. 


Assembler n’est pas un luxe réservé aux 1BB, cela 
marche très bien sur 1CC. Toutefois, l’adresse du 
programme vérifiant s’il y a un élément dans la pile 
est différente, et le programme -LEX a un code 
différent. Sur version 1CC, cette adresse est #1C4A1, 
qu’il faut donc corriger dans -+LEx (CF JPC 51 page 
23). 


Pour entrer des programmes en langage machine sur 
HP-28, la seule façon, pour le moment, est d’entrer 
au clavier les codes hexadécimaux, à coder avec ASS et 
convertir en programme avec +LEX. Toutefois, la seule 
méthode «humaine» pour réaliser une telle 
opération est l'emploi des sommes de contrôle, 
comme on le fait, quelques pages plus loin, dans le 
coin des Lhex. C’est pourquoi je préconise l’utilisation 
d’un nouveau programme ass, dont la structure 
principale reste la même : 


« 
+ LM 
« 
HEX "1 O 1 LM SIZE 
FOR X 
#1 LM X DUP2 1 + DUP SUB 3 ROLLD 
DUP SUB + + STR+ B+R DUP ROT + 
3 ROLLD CHR + SWAP 2 
STEP R-+B #FFF AND 
» 


» 


Le programme prend comme argument une chaîne de 
codes hexadécimaux dans la pile, et retourne toujours 
la chaîne codée au niveau 2 de la pile, mais aussi la 
somme de contrôle qui est un entier binaire entre 0 et 
4095. 


Pour entrer un programme en langage machine, vous 
entrez donc la première chaîne de 16 chiffres 
hexadécimaux, vous lancez ass, vous vérifiez la 
somme de contrôle, vous faites DROP, vous entrez la 
seconde chaîne, Ass, vérifiez la somme de contrôle, 
vous faites +. La dernière chaîne de codes introduite, 
vous pouvez faire -LEX, et utiliser le programme. 


Afin d'illustrer ceci, voici une fonction que j'ai appelé 
a pour faire plus Forth, qui prend une adresse en 
binaire, et retourne le contenu binaire des cinq 
quartets suivant cette adresse : 


CON(5) #02C67 
CON(5) #1C4A1 
CON(5) #02A70 
CON(5) 6 
CON(1) © 
CON(5) #18A32 
CON(5) #02096 
REL(S) FIN 
CD1EX 

D1=C 

RO=C 

A=DATI À 

D1=A 

D1=D1+ 10 
A=DATI A 
AD1EX 


« DROP » ADR 5 + Q 


« OR » ADR 
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C=DAT1 A 
D1=A 

DAT1=C A 
C=RO 

D1=C 

A=DATO A 
DO=D0+ 5 
PC=(A) 

FIN CON(5) #02F90 


45 octets 
01234566789ABCDEF sm 


000: 76C201A4C107A206 203 
001: 0000023A8169C206 1FD 
002: 3000137135108143 0E8 
003: 1311791431331471 1B8 
004: 3114511813514216 185 
005: 4808c09F20 208 


Vous pouvez appeler cette fonction comme vous le 
voulez, personnellement, j'ai réussi à l'appeler à à 
coups de poke. 


Le programme suivant ADR, qui calcule l’adresse d’une 
fonction, utilise à : 
« LOCAL « #4FO6E à 15 + @Q 10 + à » » 


Pour déterminer l’adresse de DROP, par exemple, vous 
créez un programme commençant par DROP, par 
exemple, « DROP », et vous lancez ADR. Sur ma version 
1CC, jobtiens #17921. Si, à cette adresse, vous 
ajoutez 5, et vous exécutez à, vous obtenez, versions 
1CC, #1C4A1, l'adresse du programme vérifiant la 
présence d’un élément dans la pile. En cherchant un 
peu, vous constaterez que toutes les instructions 
demandant un élément commencent par cette 
adresse. 


Ainsi, j'invite tout ceux qui programmeront en 


langage machine à indiquer avec ADR et à comment 
trouver les bonnes adresses. 


Pour finir, voici un utilitaire pour HP-71B pour 
obtenir les chaînes de codes, ainsi que les sommes de 
contrôle, pour HP-28. 


Pour faire un programme en langage d’assemblage, 
tapez ce programme avec l'éditeur de texte du 
module Forth/Assembler précédé de : 

BIN nom du fichier objet 

CHAIN -1 


Ensuite, assemblez-le avec l’assembleur et lancez le 
programme suivant qui imprimera les codes. Pour 
ceux, qui, comme moi, ne travaillent pas dans le luxe 
que procurent les imprimantes, je conseille de faire 
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DELAY INF avant. Les codes imprimés, empressez-vous 
de détruire le fichier objet. 


Serge Vaudenay (124) 


10 DIM AS$[8] 

20 INPUT ‘Fichier BIN : !;A$ 

30 CALL CODE(A$) 

40 DESTROY A$ 

50 SUB CODE(A$) 

60 DIM B$(5],C$[16] à INTEGER A,B,C,D,L 
70  B$=ADDR$S(AS$) 

80 B$=DTHS(HTD(B$)+32) 

90 C$-PEEK$(B$,5) 

100 C$=CH[S]8CH(4,4]8C$[3,3]8CH[2,2]8C$[1,1] 
110 B$=DTH$CHTD(B$)+17) 

120  A=HTD(C$)-17 

130 PRINT A/2;! octets! 


140 PRINT 

150 PRINT ! 01234566789ABCDEF sm' 

160 PRINT 

170 FOR B=0 TO (A-1) DIV 16 

180 IF B#(A-1) DIV 16 THEN L=16 ELSE L=A-16*B 
190 C$=PEEKS(DTH$SCHTD(B$)+16*B),L) 

200 D=0 

210 FOR C=1 TO LEN(C$) DIV 2 

220 D=D+HTD(C$[2*C,2*C]&CH[2*C-1,2*C-1]) 


230 NEXT C 

240 PRINT DTH$(B)[3])&!: !; 
250 PRINT (C$&' 

260 PRINT DTH$(D)1(3] 

270 NEXT B 

280 END SUB 


ill]; 


A 


“f 


HP-41 


J.M. Baillard Algèbre linéaire 10 
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ALGEBRE LINEAIRE 


En réponse à un S.OsS. lancé dans JPC 53, je vous 
propose un programme de résolution de systèmes 
linéaires, et d’inversion de matrice sur HP-41. Ce 
programme contenant de nombreuses instructions 
synthétiques, j'en donne ensuite une version non 
synthétique, plus commode à taper. 


On peut ainsi résoudre un système de 16 équations à 
16 inconnues et même de 17 équations à 17 inconnues 
si l’on exécute le programme directement à partir du 
module X-Fonctions. 


Ce programme permettant de résoudre plusieurs 
systèmes (simultanément) ayant la même matrice 
principale, il est possible en particulier d’inverser une 
matrice carrée M d’ordre n : cela revient à résoudre nr 
systèmes linéaires ayant pour matrice principale la 
matrice M. Il faut alors introduire, après les 
coefficients de la matrice M, ceux de la matrice unité 
d'ordre n (dont les éléments sont nuls sauf la 
diagonale qui ne contient que des 1, voir l'exemple 2). 
La matrice inverse M! se trouve, en fin de procédure, 
à la place de la matrice unité. On peut ainsi inverser 
une matrice 12x12 (peut être est-il possible d’inverser 
une matrice 16x16 ou 17x17 mais ce programme ne le 
permet pas). 


Enfin, je rappelle une procédure d’introduction de 
valeurs numériques dans un tableau. 


SYSTEMES LINEAIRES 
METHODE DE GAUSS - PIVOT PARTIEL 


O1*LBL "SL" 

X=0? SF 00 STON X<>Y STO a DSE a ST+ Y DSE Y 
1,001 STOT * STO M FRC STO O ISG O * LASTX 
E2 / + + STOP SIGN STO 00 

26*LBL 14 


CLST RCL O RDN + 


31*LBL O1 


CLX RCL IND Z ABS X<=Y? GTO 02 X<>Y ENTER® + 


40*LBL 02 

ISG Z GTO 01 R° INT RCL P INT X=Y? GTO 04 
X<> L FRC ST+ Y 

52*LBL 03 

CLX RCL IND L X<> IND Y STO IND L ISGY !"" 
ISG L GTO 03 RCL 00 CHS STO 00 

64*LBL 04 

RCL IND P ST* 00 X=0? GTO 07 RCL a STO Q 


JPC 59 Page 10 


71*LBL 05 

RCL Q RCL P ST+ Y ENTER" CLX X<> IND Z 
RCL IND Ÿ / SIGN ISGY ISGZ 

83*LBL 06 

CLX RCL IND Y LASTX * ST- IND Z ISGZ "" 


ISG Y GTO 06 DSE Q GTO 05 


95*LBL 07 
RCL M ST+O ISG O INT ST/L LASTX + ST+ P 
DSE a GTO 14 RCL IND O ST* 00 RCL 00 FC? 00 
X=0? STOP 


112*LBL 08 

RCL M ENTER" INT STO O / LASTX X"2 * STO P 
RCL M RCL N DSE X " * + STO Q 

129*LBL 09 


RCL M ST- P RCL O ST+ P ST+ Q RCL IND P 
ST/ IND Q DSE O FS? 30 GTO 11 RCL IND Q SIGN 
DSE P DSE Q 


144*LBL 10 
RCL IND P LASTX * ST- 
GTO 10 GTO 09 


IND Q DSE P #‘“!"! DSE Q 


154*LBL 11 


DSE N GTO 08 RCL Q RCL M FRC + CLA BEEP END 


261 octets - SIZE : n2+cxn +1 


Codes décimaux des lignes synthétiques : 


Ligne 04 : 145,118 Ligne 116 : 145,119 
ligne 06 : 145,123 Ligne 121 : 145,120 
Ligne 07 : 151,123 Ligne 122 : 144,117 
Ligne 13 : 145,117 Ligne 123 : 144,118 
Ligne 15 : 145,119 Ligne 125 : 240 

Ligne 16 : 150,119 Ligne 128 : 145,121 
Ligne 19 : 27,18 Ligne 130 : 144,117 
Ligne 23 : 145,120 Ligne 131 : 147,120 
Ligne 28 : 144,119 Ligne 132 : 144,119 
ligne 45 : 144,120 Ligne 133 : 146,120 
Ligne 58 : 240 Ligne 134 : 146,121 
ligne 65 : 144,248 Ligne 135 : 144,248 
Ligne 69 : 144,123 Ligne 136 : 149,249 
Ligne 70 : 145,121 Ligne 137 : 151,119 
Ligne 72 : 144,121 Ligne 140 : 144,249 
Ligne 73 : 144,120 Ligne 142 : 151,120 
Ligne 90 : 240 Ligne 143 : 151,121 
Ligne 93 : 151,121 Ligne 145 : 144,248 
Ligne 96 : 144,117 Ligne 148 : 147,249 
Ligne 97 : 146,119 Ligne 149 : 151,120 
Ligne 98 : 150,119 Ligne 150 : 240 

Ligne 103: 146,120 Ligne 151 : 151,121 
Ligne 104: 151,123 Ligne 155 : 151,118 
Ligne 105: 208,0,14 Ligne 157 : 144,121 


Ligne 106: 144,244 
Ligne 113: 144,117 


Ligne 158 : 144,117 


Notes : 

- Le 610 14 de la ligne 105 est un 610 longue distance. 
- Les registres P et Q étant utilisés, ne pas interrompre 
le programme en cours d’exécution. 

- Lors de l’arrêt éventuel ligne 111, ne pas modifier le 
registre alpha. 


Mode d'emploi 


Le programme permet de calculer le déterminant 
d’une matrice, de résoudre un système linéaire de n 


équations à 7 inconnues, pourvu que la matrice 
principale du système soit la même. 


Le programme transforme le (ou les) système(s) 
linéaire(s) primitif(s) en système(s) triangulaire(s) 
supérieur(s) équivalent(s) (lignes 1 à 111). Il stoppe 
ligne 111 si le déterminant de la matrice principal est 
nul ou si si drapeau F00 est armé. Enfin, il résout le 
(ou les) système(s) en plaçant les solutions à la place 
des second membres (les coefficients de la matrice 
triangulaire restant inchangés (lignes 112 à 163). 


Les données initiales sont : n = ordre du système 
dans le registre Y. c = nombre de seconds membres 
dans le registre X. 


Les données finales sont : ddd,fff = code de contrôle 
de la solution du premier système en X et le 
déterminant du système en RO00. 


Temps d’exécution : Voici, à titre indicatif, les valeurs 
t(n;c) = temps d’exécution pour la résolution de c 
systèmes de 7 équations à 7 inconnues. 

(71) = 1 mn 205 

1(10;1) = 3 mn 155 

1(12;1) = 5 mn 105 

t(14;1) = 7 mn 505 

1(16;1) = 11 mn 105 

1(17;1) = 13 mn 155 

1(7:7) = 3 mn 20s 

(10:10) = 9 mn 105 

(12:12) = 15 mn 305 
Ces durées ont été observées sur une HP-41CX. Sur 
une HP-41CV, la vitesse est environ 7 à 8% plus 
rapide. 


Utilisation : 

- Introduire les coefficients du système, colonne après 
colonne, à partir du registre RO1. 

- Introduire 71, ENTER”. 

- Introduire c, XEQ "sr. 

Après le bip ddd,fff est la première solution. 


Pour calculer seulement le déterminant de la matrice 
carrée, enter € = Ü (pas de second membre). Dans ce 
cas, le drapeau F00 s’arme automatiquement (ligne 
03). 


Exemple 1 


Soit à résoudre le système 
Trx+5y+6z+5t= 
| 3x+2y+52z2+7t= 
[2x+ V+. Z+ Ste 
léx+4y+5z2+6t- 


5 KW rm 


On a, ici, un système de 4 équations à 4 inconnues 
donc n = 4 et c = 1, donc : SIZE 021 au moins. 
Introduire les valeurs 7, 3, 2, 6... 1, 2, 3 et 4 dans les 
registre R01 à R20. Puis 

4 ENTER” 

1 XEQ "SL, 
Après le bip, 17,020 est le code de contrôle de la 
solution. 

x = -34 se trouve donc dans le registre R17, 

y = 57 dans RI8, 

z = -16 dans R19 et 

= 10 dans R20. 


Le déterminant principal du système (qui vaut ici 1) 
se trouve dans le registre ROO. 


Exemple 2 


Trouver l'inverse de la matrice : 
[7565 
M = en 
its | 
Lé4se 


On sait que cela revient à résoudre 4 systèmes 
linéaires ayant la matrice M pour matrice principale 
et pour seconds membres : 


1]fo]fol Lo 
el lslul 
lollollal do 
LoJlollo] Î: 
On a doncicin = 4etc = 4 donc: 
- SIZE 033 au moins, 
- Introduire les 32 valeurs (colonne après colonne) 
dans les registres R01 à R32 

[75651000] 

[32570100 

[21130010 

Lé4560001| 
- Puis 4 ENTER" XEQ "SL", 
Au bip, on a le code de la première colonne de la 
matrice inverse : 17,020. 


| 
| 
| 
] 
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Les 16 coefficients de la matrice inverse M‘! se 
trouvent dans les registres R17 à R32 (à la place des 4 
seconds membres) : 
[ 19 4 21-31 
Ml l31 -7-36 51 


| 
8 2 8-13 | 
-5 -1 -5 8) 
Le déterminant de M se trouve, bien entendu, dans le 
registre ROO. 


VERSION NON SYNTHETIQUE 


O1*LBL "SL2" 

X=0? SF 00 STO 02 X<>Y STO 06 DSE 06 ST+ Y 

DSE Y 1,001 * STO 01 LASTX 6 * X<>Y FRC + 
STO 03 ISG 03 X<> L * LASTX 100 / + 7,007 + 
STO 04 SIGN STO 00 


32*LBL 16 
CLST RCL 03 RDN + 


37*LBL 01 


CLX RCL IND Z ABS X<=Y? GTO 02 X<>Y ENTER* + 


46*LBL 02 

ISG Z GTO 01 R° INT RCL 04 INT X=Y? GTO 04 
X<> L  FRC ST+ Y 

58*LBL 03 


CLX RCL IND L X<> IND Y STO IND L 
ISG L GTO 03 RCL 00 CHS STO 00 


ISG Y CLX 


7O*LBL 04 
RCL IND 04 ST* 00 X=0? GTO 07 RCL 06 STO 05 


77*LBL 05 

RCL 05 RCL 04 ST+ Y ENTER" CLX X<> IND Z 

RCL IND Y / SIGN ISG Y ISG2Z 

89*LBL 06 

CLX RCL IND Y LASTX * ST- IND Z ISG Z CLX 

ISG Y GTO 06 DSE 05 GTO 05 

101*LBL 07 

RCL 01 ST+ 03 1SG 03 INT ST/ L LASTX + ST+ 04 
DSE 06 GTO 16 RCL IND 03 ST* 00 RCL 00 FC? 00 

X=0? STOP 

118*LBL 08 

RCL 01 ENTER" INT STO 03 / LASTX X”°2 6 + * 
STO 04 RCL 01 RCL O2 1 - * + STO 05 

137*LBL 09 

RCL 01 ST- 04 RCL 03 ST+ 04 ST+ 05 RCL IND 04 


ST/ IND 05 DSE 03 FS? 30 GTO 11 RCL IND 05 SIGN 


DSE 04 DSE 05 
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152*LBL 10 
RCL IND 04 LASTX * 
DSE 05 GTO 10 GTO 09 


ST- IND 05 DSE 04 CLX 


162*LBL 11 


DSE 02 GTO 08 RCL 05 RCL 01 FRC + BEEP END 


251 octets - SIZE : n2+cxn +7 


L'utilisation de ce programme est identique à celle du 
précédent, à ceci près que les coefficients doivent être 
introduits à partir du registre R07 au lieu de RO1. 
D’autre part, le temps d’exécution est très légèrement 
plus grand. Enfin, contrairement au programme 
précédent, on peut utiliser le registre alpha (lors de 
l’éventuel arrêt ligne 117) pour visualiser la matrice 
triangulaire supérieure. 


INTRODUCTION DE VALEURS NUMERIQUES 
DANS UN TABLEAU 


Exemple 


Introduire les valeurs 6, 1/3, /2 et 7 dans un tableau à 
partir du registre R72. 


72 XEQ MTEN" + N1 = ? 
6 R/S + N2=7? 
3 1/X R/S + N3 = ? 
2 SORT R/S -+ N4 = ? 
7 R/S + N5=7? 
R/S + 72,075 (code de contrôle du tableau) 


Lorsque toutes les valeurs ont été introduites ; lors du 
PROMPT suivant, taper simplement R/S. La procédure 
prend fin en affichant le code de contrôle du tableau 
en question. 


Donnée initiale : ddd = base du tableau (ici 72) dans 
le registre X. 

Donnée finale : dddfff = le code de contrôle du 
tableau dans le registre Y (ici 72,075). 


Remarques 


-On peut remplacer les instructions RCLFLAG et 
sToFLAG des lignes 11 et 17 par les instructions 
synthétiques RCL d (code décimal : 144,126) et STo d 
(code décimal : 145,126) respectivement. 

- Si l’on veut introduire une valeur comme x dans un 
tableau, il faut armer le drapeau 22 (sans quoi la 
procédure prendrait fin) : PI, SF 22, R/S... 

- Lors des PRoMPT, la pile contient des valeurs utiles en 
Y et X qui doivent se trouver en Z et Y lors du R/S 
suivant. 


Jean-Marc Baiïllard (407) 


HP75 


E. Gengoux 


E. Gengoux 
E. Gengoux 
E. Gengoux 
E. Gengoux 
E. Gengoux 


Séries chronologiques 


Programme "LES" 
Programme "MMC" 
Programme "MMEDC" 


Programme "PREVLED" 


Programme "SPLOT" 


14 


23 
23 
23 
23 
24 
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BOITE A OUTILS POUR 
SERIES CHRONOLOGIQUES 


Les outils qui suivent sont des fonctions étendues 
VisiCalc qui ont pour but de faciliter l'étude de 
problèmes tels que : désaisonnalisation rapide d’une 
série chronologique à périodicité quelconque 
(mensuelle ou autre), estimation de sa tendance et 
prévision par extrapolation de quelques périodes 
au-delà de sa dernière valeur connue, visualisation 


graphique enfin. 


Le VisiCalc du HP-75 se révèle en ce domaine un 
instrument bien adapté, grâce à la possibilité qu’il a 
de pouvoir créer des fonctions sur mesure beaucoup 
plus puissantes que n'importe quelle macro de tableur 
du commerce, et aussi par le fait que son Basic est 
alors enrichi, et peut accéder au contenu des feuilles 
de données aussi facilement qu’à des variables. Les 
cinq programmes qui suivent n’ont certes pas la 
prétention de remplacer des logiciels plus puissants 
tels que CensusIl (estimation précise des 
composantes saisonnière, cyclique, tendancielle et 
aléatoire résiduelle d’une série) ou Box et Jenkins 
(prévision à court terme et modélisation par des 
processus stochastiques). Par contre, ils permettent 
de se faire une première idée des évolutions en cours 
et de les visualiser sous forme graphique sur une 
imprimante normale ThinkJet. En outre, l’un d’entre 
eux constitue une application directe de deux articles 
récemment parus dans JPC : celui de Daniel Atlan 
(Droite de Tuckey) et le dernier Lex de Jean-Yves 
hervé, TRILEX. 


Désaisonnalisation, moyennes et médianes 


Le premier programme, MMC, fait appel à une très 
classique moyenne mobile pour désaisonnaliser la 
série. Toutefois, pour être plus rigoureux, cette 
moyenne mobile est centrée, c’est-à-dire calculée 
pour pouvoir être affectée à la période milieu des n 
termes retenus, #7 étant la longueur du cycle 
saisonnier à filtrer. Ceci impose de prendre n-1 
termes, plus 0,5 fois chacun des termes extrêmes, 
pour ne pas « tomber à cheval » entre deux périodes. 


La moyenne mobile présente deux grands défauts : si 
elle élimine bien le facteur saisonnier, par contre elle 
contient encore le terme aléatoire (c’est-à-dire à la 
fois l'innovation (ce qui n’est pas encore apparu dans 
le schéma de décomposition antérieur - on doit donc 
le conserver et l’étudier -) et les accidents, anomalies 
et erreurs de mesure, grèves et autres effets de 
calendrier (dont on devrait pouvoir éliminer 
l'influence.) et est très sensible aux périodes 
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fortement déviantes (leur influence dure 7+1 
périodes)). Dans la pratique, ce second défaut est 
d’ailleurs le plus gênant. 


Le second programme fait donc appel, non plus à une 
moyenne mobile, mais à une médiane. Rappelons 
que, pour calculer la médiane de n nombres, il faut au 
préalable les classer, puis calculer la valeur du milieu 
(celle qui correspond à la fréquence cumulée de 0,5), 
soit le nombre médian si n est impair, et la moyenne 
des deux nombres encadrant le milieu si n est pair. 
On bénéficie alors de la propriété la plus 
remarquable de la médiane, celle d’éliminer par 
construction, sans aucun calcul supplémentaire, les 
valeurs déviantes. En effet, celles-ci restent aux 
extrémités du classement, donc influent peu sur le 
rang des observations du centre, qui servent à calculer 
la médiane. C’est la fameuse robustesse, en d’autres 
termes l’immunité au bruit. Cette fonction n’aurait 
pas pu être écrite en Basic, compte tenu du temps 
pris par le tri. Heureusement, TRILEX est arrivé, qui 
a permis de rédiger une routine presque aussi rapide 
qu’une moyenne mobile. 


Lissage et recherche de la tendance 


Les programmes LES et PREVLED font appel à une 
technique un peu différente, celle des lissages 
exponentiels. Il s’agit de moyennes mobiles un peu 
particulières, en ce sens que les observations qui y 
rentrent sont pondérées par des coefficients qui 
décroissent exponentiellement au fur et à mesure 
qu’elles sont plus éloignées (on dit que seul le passé 
récent influence vraiment la série, alors que dans une 
moyenne mobile simple, toutes les observations ont le 
même poids). Bien entendu, on limite la série aux 
termes les plus récents. 


Le lissage exponentiel simple ne peut être utilisé que 
pour des séries à tendance plate, de celles qui font le 
désespoir de votre Directeur Commercial. Mais, si 
l’on lisse à son tour le premier lissage, on s'aperçoit 
alors que lerreur systématique commise en 
remplaçant la série initiale par son LES se retrouve 
lorsqu'on remplace à son tour le premier lissage par 
le second. D’où l’idée de PREVLED, qui consiste à 
réinjecter l’erreur (LES-LED) dans LES pour en tirer 
une prévision non biaisée, ce coup-ci. 


Deux recommandations pour Vlutilisation de ces 
fonctions : tout d’abord, n’appliquez jamais 
directement un lissage sur la série brute, sauf bien 
entendu s’il s’agit d’une série de périodicité annuelle ; 
commencez par la désaisonnaliser par moyenne 
mobile (ou par médiane mobile, si vous craignez des 
anomalies difficiles à corriger autrement) ; alors 
seulement lissez cette CVS une première fois, puis le 
lissage obtenu une seconde fois, et enfin calculez la 


prévision lissée. C’est cette dernière que vous pourrez 
prolonger (au choix, par PREVLED, ou sur le 
graphique avec une règle), puis resaisonnaliser à 
laide des rapports (Série brute/PREVLED) 
correspondant à la période). Ensuite, adoptez une 
disposition de votre feuille comme ci-dessous, ou 
inspirez-vous des exemples annexés. 


A B C D E F 
S.brute K.correct. S.br.corr LES LED Prev.lissée 
nn... Here ehene memes ehe se shpessepessessssmsee 


Lancez vos lissages en prenant pour valeur de S’ à la 
lère période la valeur correspondante de X :; le 
lissage s'adaptera très rapidement (d’autant plus vite 
d'ailleurs que à (0<a<1) sera plus petit), seules les 
toutes premières périodes seront donc non 
significatives). Vous pourrez calculer, dans les autres 
colonnes (G et suivantes) vos coefficients saisonniers, 
votre série resaisonnalisée, etc. 


Visualisation graphique 


Pour préparer le terrain de la routine graphique qui 
vient à présent, voyez la disposition adoptée dans le 
second exemple. On y a éliminé l’intermédiaire de 
calcul désormais inutile qu’est le premier lissage, pour 
ne conserver que la CVS (médiane ou moyenne 
mobile) (colonne D), la tendance (lissage double) 
(colonne H), et bien entendu la série brute (colonnes 
A et C, avec possibilité de corriger la colonne A par 
coefficient en colonne B donnant C). 


SPLOT n’est plus une fonction étendue, mais un 
programme qui, dans l’environnement Basic (et non 
plus sous VisiCalc), exploite les données des colonnes 
C, D et H pour en tirer un graphe simultané des trois 
séries brutes, CVS et tendance. 


On voit mieux alors l’évolution de fond de la série, 
dont la tendance constitue un estimateur centré, à 
variance réduite. Elle est en effet nettement moins 
tourmentée, et reflète fidèlement les phases de hausse 
et de baisse, sans être prisonnière d’une forme de 
courbe plus ou moins bien adaptée, droite de 
régression ou autre. 


Le programme est issu d’un logiciel plus complet 
tournant sur 75, et comportant d’autres modules 
graphiques (HP-7470, PacScreen, Miniplotter, etc.). 
Notez que, si votre série est mensuelle et que vous 
avez mis en user labels (colonne 0) les noms des mois 
sous la forme mmm aa, il imprimera même les limites 
des années ©! Bien entendu, ïl détermine 
automatiquement son échelle. 


Conclusion et limites 


Certains aspects sortant du cadre d’un tel article ont 
été passés sous silence ; en particulier, celui du choix 
optimal de la constante de lissage « (dans les 
exemples, une valeur passe-partout de 0,2 a été 
retenue), et le problème moral posé par la formule de 
prévision, qui revient à prolonger linéairement le 
dernier arc connu de la tendance, en faisant 
l’hypothèse que les conditions locales ne changeront 
pas à lhorizon de prévision, considéré comme 
voisinage du dernier point. C’est d’ailleurs pourquoi 
la méthode est plus appropriée pour dégager une 
tendance, et que des procédés plus précises sont 
recommandés pour la prévision elle-même. 


Bibliographie 


Pour en savoir plus sur ces techniques, on pourra lire 
soit l'ouvrage de Makridakis et Wheelwright paru 
chez L.E.O. soit celui des mêmes auteurs et de V. 
Mc Gee, Forecasting Methods & Applications 
(J.Wiley, 1983), plus récent et plus complet, ainsi que 
le Que Sais-Je de Coutrot et Droesbeke (Techniques 
de Prévision), qui sont des ouvrages clairs et lisibles. 
Et sinon, rendez-vous rue de Lancry ! 


Bonnes prévisions ! 


Eric Gengoux (108) 
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AH ! VOUS ECRIVEZ 


Vous vous sentez en verve, mais vous ne savez pas sous quelle forme "l’équipe de 
rédaction" souhaite recevoir votre prose. C’est ici que se trouvent les réponses à 
vos questions. 


Dans la mesure du possible, vous devez nous envoyer vos écrits sur support 
magnétique (carte, cassette ou disquette). Soyez sans crainte, nous vous 
retournerons vos biens après copie. 


Si vous ne pouvez pas utiliser de support magnétique, ou ne pouvez vous rendre 
aux réunions, alors et alors seulement faites le sur papier. 


Que ce soit sur une feuille de papier, ou sur support magnétique, ne dépassez 
pas 50 caractères par ligne. 


Pour nous épargner du travail, insérez dans votre texte les commandes de 
formattage suivantes (et non les commandes du formatteur HP) : 

"?" centre un titre, par exemple : 
"TITRE 


"\'(CHR$(92)) marque le début et la fin d’un paragraphe. Par exemple : 


\Début de paragraphe exprimant le contenu de vos idées qui, même si vous en 
doutez, intéressera certains des membres du Club. Surtout si vous vous sentez 
débutant. Les articles pour débutants écrits par des débutants sont ceux qui 
manquent le plus. Fin de paragraphe.\ 


N'oubliez pas de mettre les accents. Utilisez le jeu de caractères Romans. Les 


possesseurs de HP71 utiliseront les redéfinitions de touches ci-dessous, ainsi que 
le fichier CHARLEX listé dans le coin des Lhex. 


Jean-Jacques Dhénin (177) 


DEF KEY ‘fW'!, CHR$(197); (é) 
DEF KEY 'fE', CHR$(193); (é) 
DEF KEY ‘fR', CHR$(201); (è) 
DEF KEY ‘fY!, CHR$S(203); çà) 
DEF KEY ‘fU', CHR$(195); çà) 
DEF KEY 'fl', CHR$(209); [&b) 
DEF KEY 'fO', CHR$(194); (ô) 
DEF KEY 'f/!, CHR$(92); (V0 
DEF KEY ‘fA!, CHR$(192); (à) 
DEF KEY 'fS', CHR$(200); (à) 
DEF KEY 'fD', CHR$(205); (ë) 
DEF KEY 'fJ', CHR$(207); çü) 
DEF KEY ‘fK!, CHR$S(221); ci) 
DEF KEY ‘f*1, CHR$(124); Cl) 


DEF KEY 'fC', CHR$(181); (ç) 


PPC PARIS SE REUNIT 
UNE FOIS PAR MOIS 


Comme vous le savez peut être déjà, PPC Paris se réunit une fois par mois, en 
plein coeur de Paris. Amenez votre matériel, votre bonne volonté et vos idées ! 
Plus vous en apporterez, et plus vous en trouverez chez vos collègues de PPC. 


Ces réunions se déroulent de manière très libre, aucun ordre du jour, discussion 
ou autre n'étant imposé. Un membre du bureau est toujours présent. Ainsi, si 
vous désirez remettre votre article tout frais au Journal, si vous avez des 
suggestions à faire, si vous voulez vous procurer des anciens numéros de JPC, ce 
sera en principe toujours possible. 


Si donc cela vous intéresse, n’hésitez plus un seul instant, venez nous rejoindre 
tous les premiers samedis de chaque mois (sauf en période de vacances 
scolaires) au : 


Centre de Jeunesse et de Loisirs Jean Verdier 
11 rue de Lancry 
75010 Paris 


et en montant au deuxième étage, vous entendrez des éclats de rire et des 
discussions passionnées vers la salle 215. Attention, toutefois, de venir entre 16 et 
19h. 


Pour l’accès en métro, trois possibilités s’offrent à vous : 
- Métro Strasbourg Saint Denis : 

Sortie porte St Martin / Bd St Denis, coté pairs 

- Métro République : 

Sortie Bd St Martin, coté pairs 

- Métro Jacques Bonsergent : 

Sortie Bd Magenta, coté impairs. 


Ab, j'oubliais ! JPC est (souvent) distribué en avant première lors de ces 
réunions... À bon entendeur, salut ! 


Les dates des prochaines réunions sont : 
Samedi 5 novembre 1988 

Samedi 3 décembre 1988 

Samedi 7 janvier 1989 

Samedi 4 février 1989 

Samedi 4 mars 1989 

Samedi 6 mai 1989 

Samedi 3 juin 1989 


Pierre David (37) 


eme 


NOUS EN AVONS 


La coopérative du Club vous propose : 


- de lecteurs de cartes magnétiques pour HP-71, neufs, dans leur boîte d’origine, avec 5 cartes magnétiques, 
pour 500 F (port compris), 

- des anciens numéros de JPC, au prix de 40 F + 7,40 F de frais d’affranchissement, 

- d’une année complète de numéros de JPC (février à janvier) pour 300 F (offre spéciale) port compris, 

- des I.D.S. du module Forth / Assembleur (listing interne commenté par HP) pour 250 F (port compris), 

- des VASM pour HP-41 (listings des Roms internes commenté par HP) pour 300 F (port compris), 

- de manuels de service du HP-41 au prix de 75 F (port compris), 

- de manuels de service du HP-75 au prix de 75 F (port compris). 


En outre, le module JPC Rom pour HP-71 est disponible. Vous nous adressez votre Eprom CMT (de 
préférence 64 Ko), et nous la programmons suivant une des options ci-dessous : 

- JPC Rom + Manuel, pour 600 F, 

- JPC Rom + Manuel + vos propres programmes, pour 800 F. 


Si vous souhaitez des renseignements complémentaires, n’hésitez pas à nous contacter. 


VOUS EN VOULEZ 
Nom : 
Prénom : 
No de membre : 
Adresse : 
Commande : 

Qté Prix Unitaire Prix Total 

A A ART Re AT 
anciens numéros de JPC x 47,40 FF 
année complète de JPC x 300 FF 
1.D.S. du module Forth x 250 FF 
Programmathèque HP-71 (joindre 3 disquettes) x 75 FF 
Manuel de service pour HP-41 x 75 FF 
Manuel de service pour HP-75 x 75 FF 
JPC Rom + Manuel x 600 FF 
JPC Rom + Manuel + vos propres programmes x 800 FF 
Actualisation Eprom x 150 FF 
EN A EE APR Rs tie 1e me A A FE IAE 


Préciser éventuellement Les 
numéros de JPC commandés : 


PPC PARIS 


Association régie par la loi de 1901, enregistrée 
à Paris le 2 décembre 1982 sous le numéro 82/3240 


BULLETIN D’ADHESION 
Nour EE RAR RER RO SE TE RL PO A TS ECO 
Prénom |_[_|_|_1_1_[_1_1_1_1_1_1_1_1_1_11 Date de naissance |_|_1/1_1_1/1_1_] 
a rh BU PA CE ER NA OP EL RE A EP A, 
POPÉRORE PETER SEINE LORIE OO SR OT TEEN 
ADI RP RE ERA AO ESP SRIOETL ESS 
LE PAS PTE EAP A PAP OA ee 
Gode Postal: | ILILELE MIE LILAS PE MIRE ER RREeEEOr CORNE 
1-1 EAP PE LR EURO 


Profession 
Intérêts 


Matériel HP en votre possession 


Autre matériel informatique 


Comment avez-vous connu PPC Paris ? 


Que recherchez-vous au sein de PPC Paris ? 


La Loi No 78-17 du 6 janvier 1978 relative à l'informatique, aux fichiers et 
aux Libertés, garantit à toute personne justifiant de son identité un droit 
d'accès et de rectification auprès des services ou organismes chargés de mettre 
en oeuvre des traitements informatiques comportant des informations nominatives 
Le concernant. 


Je souhaite adhérer au club PPC Paris conformément aux statuts de 
l'Association. AU mieux de ma connaissance, je déclare avoir Le droit de 
fournir tous Les programmes et informations que je vous enverrai (sans 
enfreindre des obligations de secret à l'égard d'autres personnes ou 
organismes) pour publication dans Le Journal de Liaison, sans obligations ni 
responsabilité d'aucune sorte (en cas d'utilisation frauduleuse) de La part des 
dirigeants de PPC-Paris. 


Date |_|_[/1_1_1/191 1] 
Signature, précédée de La mention "Lu et Approuvé!! 


Le montant de La cotisation s'élève à 350.00 F pour un an. 
Etudiants: 300.00 F (justificatif indispensable) 


Paiement à l'ordre de "PPC Paris! par chèque bancaire ou virement postal (CCP 
No 18 823 40 C à Paris). Ne pas utiliser d'Eurochèques. 


Veuillez envoyer toute correspondance à : 
PPC Paris, BP 604, 75028 Paris Cedex 01, France 


ASSEMBLEUR 


G. Toublanc 


BASIC 
B. Gil 


E. Gengoux 
E. Gengoux 
E. Gengoux 
E. Gengoux 
E. Gengoux 
G. Toublanc 
B. Gil 


LE COIN DES LHEX 


Graphes de Connett (acte IT) 


Tracé de courbes 


Programme "LES" pour HP-75 
Programme "MMC" pour HP-75 
Programme "MMEDC" pour HP-75 
Programme "PREVLED" pour HP-75 
Programme "SPLOT" pour HP-75 
Programme "CONNETT" 
Programme "COURBE!" 
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GRAPHES DE CONNETT 
(ACTE Il) 


Stefano Tendon nous offrait dans JPC44 un 
programme Basic réalisant des graphes de Connett 
avec une imprimante ThinkJet. Une optimisation 
avait permis à l’auteur un gain de rapidité de 225 %. 


Malgré ce bon résultat il fallait attendre 6 heures et 
40 minutes avant d'admirer le graphe exemple de la 
page 29 de JPC 44. 


En remplaçant la boucle interne du programme Basic 
par une routine en assembleur il est possible 
d'obtenir le résultat 15,7 fois plus vite, c’est à dire en 
25 minutes et demie au lieu de 400 minutes. 


En prime, je vous offre une deuxième fonction 
obtenue en compliquant la formule mais donnant des 
graphes assez inattendus selon les paramètres choisis. 


Paramètres à préciser 


X,Y : coordonnées exprimées par 2 réels 

R : résolution 0<R < 640 

L : nombre d'unités O0 <L<R 

Toutes ces données sont converties en nombres 
entiers si elles ne le sont pas déjà. Après arrondi 
R<>0,. 


Option : 
soit graphe de Connett + C 
soit graphe dérivé + D 
dans ce cas X<0 et Y>0 


Il ne vous reste plus qu’à découvrir les bonnes 
combinaisons de paramètres donnant les plus beaux 
graphes. 


Guy Toublanc (276) 
LEX !CONNLEX! 
t EQU 145 
ID #5C 
MSG 0 
POLL O0 


ENTRY  CONSe 
CHAR  #D 
ENTRY  GRAPe 
CHAR  #D 

KEY !CONPS'! 
TOKEN t 
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ADRS50 
ARGERR 
COMCK 
DCHXW 
DROPDC 
EXPEX- 
IDIV 
MFERR 
MPY 
MSPARe 
NUMCK 
NXTSTM 
POPIÎR 
RESPTR 
RJUST 
RNDAHX 


popir 


L37 


popa 


argerr 


CONSd 
CONSp 
L40 


resptr 
mspare 
invvar 


GRAPe 


KEY 
TOKEN 
ENDTXT 


EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 
EQU 


D1=D1+ 
GOSBVL 
ST=0 
?A=0 
GOYES 
ST=1 
GOSBVL 
C=A 
GOVLNG 
D1=D1+ 
GOSBVL 
GONC 
c=0 
C=A 
RTN 
GOVLNG 
GOVLNG 
LC(1) 
C=Cc-1 
R2=C 
GOSBVL 
GOSUB 
C=R2 
?2C=0 
GOYES 
GOSBVL 
GONC 
C=R2 
GOC 
GOVLNG 
GOVLNG 
LC(2) 
GOVLNG 
REL(5) 
REL(5) 
ST=1 
GOTO 


'GRAPS! 
dt 


#0F551 
#0BF19 
#036CD 
#OECDC 
#05470 
#0F178 
#O0EC7B 
#09393 
#OECBB 
#O2E5C 
#0369D 
#08A48 
#OE8FD 
#03172 
#12AE2 
#136CB 


16 
POPIR 
0 

S 

L37 

0 
RJUST 
W 
DCHXW 
16 
RNDAHX 
argerr 
W 

A 


ARGERR 
DROPDC 
6 
p 


NUMCK 
resptr 


p 
resptr 
COMCK 
mspare 


L40 
RESPTR 
MSPARe 
83 
MFERR 
CONSd 
CONSp 
3 

CONS 


CONSe 
CONS 


BOUCL 


L94 


L99 


ADD 
MUL 


REL(5) 
REL(5) 
ST=0 
GOSBVL 
LC(3) 
B=C 
GOSBVL 
Goc 
DO=D0+ 
C=DATO 
ADOEX 
A=A-C 
DO=A 
DO=D0+ 
LC(1) 
B=C 
c=0 
DATO=C 
DO=D0+ 
B=B-1 
GONC 
GOSBVL 
ST=0 
?A=0 
GOYES 
ST=1 
GOSUB 
ST=0 
?2ST=0 
GOYES 
ST=1 
RO=C 
GOSUB 
R1=C 
GOSUB 
R2=C 
GOSUB 
R3=C 
A=R2 
GOSBVL 
A=RO 
?2ST=0 
GOYES 
A=A-C 
GONC 
A=-A 
ST=0 
GOC 
A=A+C 
C=R3 
GOSBVL 
RO=A 
GOSUB 
D=C 
A=C 
GOSBVL 
R3=A 
A=R2 


CONSd 
CONSp 
3 
EXPEX- 
‘g! 

X 
ADRS50 
invvar 
11 

A 


> 


ADD2 
STO 


BOUCL 
POPIR 
2 

M 

L94 


popir 


L99 


Ps 
popir PS2 
popa 


popa 


MPY 


ADD 
W L138 
MUL L139 


MUL 


MPY 


popa 
W 
W L149 
MPY 
L151 


C=D 
C=C+1 
GOSBVL 
A=R1 
2ST=0 
GOYES 
A=A-C 
GONC 
A=-A 
ST=0 
GOC 
A=A+C 
R1=A 
C=D 
A=C 
A=A-1 
B=A 
ASRB 
ASRB 
LC(1) 
B=B&C 
B=C-B 
?B=0 
GOYES 
A=A+1 
CDOEX 
C=C-A 
DO=C 
D=0 
D=D+1 
GONC 
D=D+D 
B=B-1 
GONC 
A=R1 
C=R2 
?2ST=1 
GOYES 
A=A-C 
GONC 
A=-A 
ST=1 
GOC 
A=A+C 
R1=A 
C=D 
GOSBVL 
A=RO 
2ST=1 
GOYES 
?2ST=0 
GOYES 
?ST=0 
GOYES 
?ST=1 
GOYES 
A=A-C 
GONC 


W 
A 
MPY 


ADD2 


STO 


STO 


V»>£ZZ 


L138 


L139 


L139 


MPY 


L149 
L155 
L151 
L155 


L156 
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L155 
L156 


L162 


DCR 


DCR2 


FIN 
PSET 


A=-A VW 
GOC L156 
A=A+C VW 
C=R3 

GOSBVL IDIV 
?ST=0 2 
GOYES L162 
A=A+1 B 
SB=0 

ASRB 

?SB=0 

GOYES PSET 
D=D+D S 
GONC  DCR2 
D=D+1 S 
DO=DO+ 1 
D=D-1 X 
?2D#0 X 
GOYES  LOOP 
GOVLNG NXTSTM 
C=DATO S 
C=CID S 
DATO=C S 
GOC DCR 


END 
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TZ à ZA. 


SS 


\ 
\ 


4, 
Î 
' 


RSS 
ETS, 
=2- 2 


MS ==. è x 





D 631 -9091 82646261 49 
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TRACE DE COURBES 


Le programme que je vous propose permet de tracer 
des courbes sur le HP-71. Ce problème a déjà été 
traité dans JPC, mais à la différence du précédent 
programme qui utilisait le module Forth/Assembleur, 
celui-ci est entièrement en Basic. Il utilise toutefois le 
très utile Lex STRUC2 et le module Maths. 


Mode d’emploi 


Vous devez d’abord indiquer s’il s’agit d’une courbe 
de type x(f), y(?) (choix 1) ou y = f(x) (choix 2). 


Le premier choix permet de tracer des courbes en 
paramètré et en polaire (r = f(t) donne x(t) = f(t)cos t 
et y(f) = f(tsin 1) ce qui donnera la même courbe. 
Dans cette partie le repère est orthonormé. Le 
nombre de points permet d’ajuster la précision du 
tracé désirée. S’il est trop grand prenez votre mal en 
patience. 





Le deuxième choix trace des fonctions mais le repère 
n’est pas orthonormé : il utilise au maximum l'écran 
en calculant l'échelle. 


L'écran est de 64x132 points. L’augmenter n’est pas 
un problème (voir listing) : changer les échelles et 
ajouter des A$C..) [132]. 


Pour visualiser la courbe, utilisez les touches de 
curseur [*] et [v]. La touche [ON] permet d’arrêter et 
[ENDLINE] revient au menu. 


L’inconvénient majeur de mon programme est que le 


déplacement est de 8 lignes en 8 lignes. Dès que j'en 
aurai le temps, j’essaierai d’en faire un qui se déplace 


de 4 en 4. 
Bruno Gil (332) cb 
o 





Dretkes | 
aauches | 









Git de demarrage entrainant un octet poursuivit par een 
bit de parité … 
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Programme LES (pour HP-75 et Rom VisiCalc) 


- LES(<alpha>,S'(t-1),X(t-1)) Liss.expon.simple 


20 
30 
40 
50 
60 
70 
80 


ASSIGN # 255 TO "VISIDATA" 
READ # 255,0 ; N1,N2 

IF N1#3 OR N2#3 THEN END 
READ # 255 ; A,S1,X 

IF A<O OR A>1 THEN END 
S=A*X+(1-A)*S1 

PRINT # 255,0 ; S 


he ee he he ee ee he ee eee ee he he ee eee he eh eee eh ee ee he he eee oh eh eh ee he eh eh he ee ee ee RER RAR 


Programme MMC (pour HP-75 et Rom VisiCalc) 


20 
30 
40 
50 
60 
70 
80 
90 
100 
110 


MMC(Liste) (max 0.5+11+0.5 périodes) 
Moyenne mobile centrée 
ASSIGN # 255 TO "VISIDATA'" 
READ # 255,0 ; N1,N2 
IF FPCN1/2)=0 THEN END 
IF N1#N2 OR N1>13 THEN END 
READ # 255 ; M à M=M/2 
FOR 1=2 TO N1-1 

READ # 255 ; M1 @ M=M+M1 
NEXT 1! 
READ # 255 ; M1 à M=M+M1/2 
M=M/(N1-1) @ PRINT # 255,0 ; M 


ee he he he he eh eh he eh he he he eh he ee he eh ee ee he he eh A ee RAA ARR 


Programme MMEDC (pour HP-75 et Rom VisiCalc) 


20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 


MMEDC(lListe) médiane mobile centrée 
ASSIGN # 255 TO "VISIDATA" 

OPTION BASE 1 

READ # 255,0 ; N1,N2 

IF N1>13 OR N1#N2 THEN END 

DIM T(13) 

FOR 1=1 TO N1 

READ # 255 ; T(I1) 

NEXT I 

TRITAB T(1),1 

IF MOD(N1,2)#0 THEN M=T(N1/2+.5) à GOTO 130 
M=.5*(T(N1/2)+TC1+N1/2)) 

PRINT # 255,0 ; M 


VOHONONO HONOR HOMO AURA AR Re he RO AR 


Programme PREVLED (pour HP-75 et Rom VisiCalc) 
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- PREVLED(<horizon>,<alpha>,S'(t},S"(t)) 

Prévision par Liss.expon.double 

20 ASSIGN # 255 TO "VISIDATA!" 

30 INTEGER M 

40 READ # 255,0 ; N1,N2 

50 IF N1#4 OR N2#4 THEN END 

60 READ # 255 ; M,A,S1,S2 

70 IF M<1 OR A<0 OR A>1 THEN END 

80 A1=2*S1-S2 à B1=A*(S1-S2)/(1-A) 

90 S=A1+M*B1 

100 PRINT # 255,0 ; S 


AR AEHHHEH RENE Ke HER HD HR AH HA HA A A A A AO AR OH AN EEE ee 


Programme SPLOT (pour HP-75 et Rom VisiCalc) 


- SPLOT [imprimante Ligne] (c)EG V3.0 
20 OPTION BASE 1 
30 DIM W$[8],M$[8],G$[132] 
40 PRINT CHR$(27)&'&k2S"! à PWIDTH 132 
50 PRINT ‘Graphe série sous VC! à PRINT à PRINT 
60 INPUT "Fichier VC : "; W$ à W$=UPRCS(WS) 
70 IF INCAT(W$,"W!")#2 THEN BEEP à GOTO 60 
80 WORKSHEET W$ à M=1 
90 INPUT f'#ligne 1ere periode? ",111:M0 
100 IF MO<1 OR MO>MAXROW THEN BEEP à GOTO 90 
110 INPUT "#ligne derniere per.?",STR$(MAXROW) ; M9 
120 IF M9>MAXROW OR M9<M1 THEN BEEP @ GOTO 110 
130 INPUT "Nbre unites/module=?","100";14 
140 IF MOD(14,100)#0 THEN BEEP à GOTO 130 
150 DISP "Tracer Le GRID O/N ?" 
160 K$=UPRCS(WAITKEY$(-1)) à M2=POS(!NO!,K$) 
170 IF NOT M2 THEN BEEP à GOTO 160 ELSE M2=M2-1 
180 DISP "Un instant..." 
190 BO=VAL(GETVALUES(3,MAXROW\2,1)) à B9=B0 
200 AO=VAL(GETVALUES(4,MAXROW\2,1)) @ A9=A0 
210 CO=VAL(GETVALUES(5,MAXROW\2,1)) à C9=Cc0 
220 FOR 1=1 TO MAXROW ! Recherche mini/maxi 
230 A=VAL(GETVALUES(4,1,1)) 
240 B=VAL(GETVALUES(3,1,1)) 
250 C=VAL(GETVALUES(5,1,1)) 
260 IF A<AO AND A#O THEN AO=A 
270 IF A>A9 AND A#9999 THEN A9=A 
280 IF B<BO AND B#0 THEN B0=B 
290 IF B>B9 AND B#9999 THEN B9=B 


300 IF C<CO AND C#0 THEN CO=C 
310 IF C>C9 AND C#9999 THEN C9=C 
320 NEXT 1! 


330 AO=MIN(AO,B0) à A9=MAX(A9,B9) ! Choix échelle 

340 I=CEIL((A9-A0)/14) à 12=CEIL(100/1) 

350 PRINT "Désaisonnalisation de La série : ",4$ 

360 PRINT “Minimum= ";B0," Maximum= ";B9 

370 PRINT "'(s) série brute, (+) série MMobile, (*) tendance." à PRINT 
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380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 


FOR K=0 TO 1 
PRINT TAB(K*12+4);1P(10*(A0+14*K))/10; 
NEXT K 
GOSUB 640 à GOSUB 660 à PRINT G$ 
PRINT 


M$=GETLABELS$(0,1) à IF M$[1,3]="DEC" THEN GOTO 450 


GOSUB 640 à GOSUB 660 à PRINT G$ 
FOR 11=M0 TO M9 
MS$=GETLABELS(O,11)8& 1" 
IF M$[1,3]="JAN" THEN GOSUB 640 
IF M2=1 THEN GOSUB 660 
G$[1,81=M$ 
A1=VAL(GETVALUES(3,11,1)) 
IF A1=0 OR A1=9999 THEN GOTO 530 
A1=IPCCAT-A0)*12/14)+8 à GSCA1,A1]="s" 
A2=VAL(GETVALUES(4,11,1)) 
IF M=1 THEN GOTO 570 
IF A2=0 OR A2=9999 THEN GOTO 570 
A2=IPCCA2-A0)*12/14)+8 à GS$[A2,A2]="+1 
A3=VAL(GETVALUES(7,11,1)) 
IF A3=0 OR A3=9999 THEN GOTO 600 
A3=IPCCA3-A0)*12/14)+8 à GS[A3,A3]="*1 
PRINT G$ à GOSUB 670 à NEXT 11 
MS$=GETLABEL$(0,MAXROW) 


IF M$[1,3]#"JAN" THEN GOSUB 640 à GOSUB 660 à PRINT G$ 


BEEP à DISP "Termine!' Q END 


G$[1,8]=" " 

FOR K=9 TO I*12+8 à G$[K,K]="-1" à NEXT K @ RETURN 

FOR K=0 TO I @ G$[8+K*12,8+K*12]=1"I" @ NEXT K @ RETURN 
FOR K=9 TO I*12+8 àQ G$[K,K]=" " @ NEXT K @ RETURN 


Ve he che he ee he he de he eee ee ee eh ee eh He A AA ee ee AAA RAA OR ARR 


Programme "CONNETT" (tracé de graphes sur imprimante ThinkJet. Nécessite CONNLEX) 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 


140 


Graphes de Connett avec l'imprimante Thinkjet 

d'après le programme basic de Stéphano Tendon (JPC 44) 
nécessite Le Lex CONNLEX pour Les mots CONPS et GRAPS 
modifications du programmme basic et conception 


du Lex CONNLEX : auteur Guy Toublanc (276) 


PRINTER IS :PRINTER @ PWIDTH INF 
DIM X,Y,L,A,B,C,D,B,1,R,G$1(80] ,R$[1] 
G$=CHRS(O)ECHRE(O)ECHRE(O )ECHRE(O)ECHRE(O) 
FOR 1=1 TO 4 

G$=G$&GS 
NEXT I 
INPUT MR, X, Y, L 2? MR,X,Y,L 
INPUT ‘“Connett/Derive C/D ? ";:R$ 
PRINT CHR$(27)&"*r640sA!; 
A=MOD (X*X+Y*Y,2) 
B=2*R*L*X @ C=2*R*L*Y @ D=L*L 
IF R$="D'" THEN 180 
Boucle pour l'option graphe de Connett 
FOR 1=1 TOR 
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150 CONPS R,1,D,C,B,A @ PRINT CHR$(27)&"*bB80W"8G$; 
160 NEXT I à GOTO 210 
- Boucle pour l'option graphe dérivé du graphe de Connett 
180 FOR I=1 TOR 
190 GRAPS R,1,D,C,B,A @ PRINT CHR$(27)&"*bD80W"&GS; 
200 NEXT 1! 
210 PRINT CHR$(27)&"*rB! 
220 PRINT R;X;Y;L @ PRINT 


AOHH HA He He He He He He HR DH HR He A HA HA AR AAA AA OH ee ee ee 


Programme "COURBET" (tracé de courbes sur l'écran. Nécessite STRUC2 et module Maths) 


10 DESTROY ALL 

20 DIM A$(8) [132] ,F$[132] 
30 FOR 1=1 TO 132 

40 F$=F$&CHRS(0) 

50 NEXT ! 

60 FOR 1=1 TO 8 

70 ASC1)=F$ 

80 NEXT ! 

90 A=0 
100 INPUT 11) x(t),yCt), 2) y=f(x) ? !;:B$ 
110 IF VAL(B$)=1 THEN ‘PARA! 
120 IF VAL(B$)=2 THEN 'CART' 
130 BEEP à GOTO 100 
140 IF Y<=8 AND Y>0 THEN 


150 A=1 

160 GOSUB 'CONVERT! 
170 END IF 

180 IF Y<=16 AND Y>8 THEN 
190 A=2 

200 Y=Y-8 

210 GOSUB 'CONVERT! 
220 END IF 

230 IF Y<=24 AND Y>16 THEN 
240 A=3 

250 Y=Y-16 

260 GOSUB '"CONVERT'! 
270 END IF 

280 IF Y<=32 AND Y>24 THEN 
290 A=4 


300 Y=Y-24 
310 GOSUB !CONVERT! 


320 END IF 

330 IF Y<=40 AND Y>32 THEN 
340 A=5 

350 Y=Y-32 

360 GOSUB 'CONVERT'! 
370 END IF 

380 IF Y<=48 AND Y>40 THEN 
390 A=6 

400 Y=Y-40 

410 GOSUB ‘'CONVERT'! 
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420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 


540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 


790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 


910 
920 
930 
940 
950 
960 


END IF 
IF Y<=56 AND Y>48 THEN 
A=7 
Y=Y-48 
GOSUB !CONVERT! 
END IF 
IF Y<=64 AND Y>56 THEN 
A=8 
Y=Y-56 
GOSUB !CONVERT' 
END IF 
RETURN 


"AFFICH!': 
BEEP 
DISP ‘Affichage de la courbe! 
T=1 
GDISP A$(T) 
IF T=4 THEN SFLAG 0 ELSE CFLAG 0 
REPEAT 
K$S=KEYS 
UNTIL K$#!! 
SELECT K$ 
CASE ‘#50! 
T=T#1 
IF T>=8 THEN T=8 
CASE ‘#51! 
T=7-1 
IF T<=1 THEN T=1 
CASE ‘#38! 
R=1 
CASE ‘#43! 
STOP 
END SELECT 
IF R=1 THEN R=0 àQ GOTO 60 
GOTO 580 
CFLAG 0 
STOP 


*CONVERT  : 
Z=NUMCASCA)[T, T1) 
IF Z=0 THEN 
S$= 00000000! 
ELSE 
S$=BSTRS(Z,2) 
FOR C=1 TO 8-LEN(S$) à S$='0'8S$ à NEXT C 
END IF 
S$IY,YI="1" 
Z=BVAL(S$,2) 


ASCA)LT,TI=CHR$(2) 
RETURN 

"CART!': 

E=0 

INPUT 'F(X)=';HS$ 


INPUT 'Xmini !;X1 
INPUT 'Xmaxi !;X2 
M=0 
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970 

980 

990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 


1140 
1150 
1160 
1170 
1180 
1190 
1200 


1210 
1220 
1230 
1240 
1250 


1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 


IF X1>X2 THEN !CART!' 

FOR T=1 TO 132 
X=X1+(X2-X1)/131*(T-1) 
Y=VALCH$) @ Y=INT(Y) 
IF ABS(Y)>M THEN M=ABS(Y) 

NEXT T 

FOR T=1 TO 132 
X=X1+(X2-X1)/131*CT-1) 
IF ABS(X)<.05 THEN E=T 
Y=VALCHS) 
Y=INT(Y*31/M)+32 
GOSUB 140 

NEXT T 

GOSUB 'AXEY' 

IF E>0 THEN GOSUB 'AXEX! 

GOTO 'AFFICH! 

STOP 


"AXEY!: 

A=4 à Y=8 

FOR T=1 TO 132 
GOSUB !CONVERT! 

NEXT T 

RETURN 

STOP 


'AXEX!: 

FOR 1=1 TO 8 
ASCI)LE,EJ=CHR$S(255) 

NEXT 1! 

RETURN 


"PARA! : 
E=0 
INPUT 'XCT)= !;:X$ 
INPUT 'YCT)= 5: Y$ 
INPUT 'tmini !;T1 
INPUT 'tmaxi !;T2 
INPUT Nb points 2!,!100!;V 
IF T1>T2 THEN 'CART!' 
M1=0 à M2=0 à M=0 
FOR T=T1 TO T2 STEP (T2-T1)/V 
X=VAL(XS) 
Y=VAL(Y$) 
IF ABS(X)>M1 THEN M1=ABS(X) 
IF ABS(Y)>M2 THEN M2=ABS(Y) 
NEXT T 
IF 65/M1>32/M2 THEN 
M=31/M2 
ELSE 
M=65/M1 
END IF 
FOR U=T1 TO T2 STEP (T2-T1)/V 
T=U 
X=INT(VAL(X$)*M)+66 
Y=INT(VAL(YS)*M)+32 
T=X 
GOSUB 140 
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1520 NEXT U 

1530 GOSUB 'AXEY' 

1540 E=66 à GOSUB 'AXEX' 
1550 GOTO 'AFFICH! 
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LE COIN DES LHEX 


Comme de coutume, cette rubrique contient la liste des codes hexadécimaux des fichiers Lex 
parus ce mois-ci. 


Rappelons ce qu’est un fichier Lex : c’est un programme pour le HP-71, en assembleur, qui 
apporte de nouvelles fonctions. Celles-ci sont utilisables directement, ou dans des programmes 
Basic. 


Pour bénéficier de ces nouvelles fonctions, vous n’avez pas besoin de programmer vous-même en 
assembleur, ni de posséder un module Forth/Assembleur. 


Il suffit de recopier le petit programme basic "MAKELEX" ci-dessous, de le lancer et de recopier 
les codes du fichier Lex désiré. Quand vous avez fini, les nouvelles fonctions sont accessibles, 
après avoir éteint et rallumé votre HP-71. 


Si l'erreur "Erreur de somme" apparaît, vérifiez la ligne que vous avez introduite. 
Vous trouverez donc le Lex CHARLEX nécessaire à la rédaction de votre article (voir "Ah ! Vous 


écrivez !"), ainsi que le Lex CONNLEX de Guy Toublanc. En outre, le programme COURBE! de 
Bruno Gil nécessite le Lex de programmation structurée STRUC2. 


CHARLEX 

CONNLEX  CONPS XWORD 92145 GRAPS XWORD 92146 

STRUC2 END XWORD 225066 WHILE XWORD 225067 
REPEAT XWORD 225068 UNTIL XWORD 225069 
LEAVE XWORD 225070 LOOP XWORD 225096 
SELECT XWORD 225097 CASE XWORD 225098 
IF XWORD 225099 ELSE XWORD 225100 


10 CALL MLEX à SUB MLEX à SFLAG -1 à PURGE AH à INPUT "Nb. d'octets: ":N @ LC OFF 

20 CREATE DATA AH,1,N-4 @ A=HTD(ADDRS("AH")) à B=A @ GOSUB 130 

30 Q=1 à X=0 à INPUT "000: ",P$;A$ à C$=A$ à S=0 à GOSUB 90 

40 Q=2 à X=1 à GOSUB 80 à AS=A$8CS à A=A+37 à N=N*2+37 à Q=3 à SFLAG 5 à FOR X=2 TO N DIV 16-1 
50 GOSUB 80 à C$=C$[5*FLAG(5)+1] à POKE DTHS(A),CS @ A=A+16-5*FLAG(5,0) à NEXT X à Q=4 
60 DISP DTH$CX)(3]; à INPUT ": ",P$[1,MOD(N,16)];C$ à GOSUB 90 

70 POKE DTHS(A),C$ à POKE DTH$(B),A$ à CFLAG -1 à END 

80 DISP DTH$CX)[3]; @ INPUT ": ,P$;C$ 

90 DISP DTH$CX)[3]; à INPUT " sm ","---11:D$ 

100 M=S à FOR Z=1 TO LEN(C$) à M=NUM(C$[Z])+M+1 @ NEXT Z 

110 IF D$=DTH$(MOD(M,4096))13] THEN GOSUB 130 à S=M à RETURN 

120 DISP "Erreur de somme" à BEEP à P$-C$ à POP 4 ON Q GOTO 30,40,50,60 

130 P$=U------.-........ 1 à RETURN 
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CHARLEX 


000: 
001: 
002: 
003: 
004: 
005: 
006: 
007: 
008: 
009: 
00A: 
00B: 
00C: 
00D: 
00E : 
00F: 
010: 
011: 
012: 
013: 
014: 
015: 
016: 
017: 
018: 
019: 
01A: 
018: 
01C: 
010: 
01E: 
01F: 
020: 
021: 
022: 
023: 
024: 
025: 
026: 
027: 
028: 
029: 
O2A: 
028: 
02c: 
02: 
O2E: 
02F: 
030: 
031: 
032: 
033: 
034: 
035: 


0123456789ABCDEF 


34841425C4548502 
802E000000000000 
5E4001EFF0000000 
FE0000000800001F 
F31BF961400032BF 
38F14A11DB10AD23 
07D532BFB8FD7911 
11AD754D7A101743 
11014D1CB15D0000 
71450375FF864834 
5655581008355654 
5810070507701724 
7700775070077517 
2077040708364545 
4A30000449724000 
0808094A2C180814 
A464242008355455 
581000054c714000 
0C3142404C700832 
41414A70002078A0 
2F30000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
000000000000080€ 
1A28080008080A2C 
180008040€ 340800 
08001E3018000000 
0000000000000000 
0000000000000000 
0000000000000000 
0201000000010200 
0000000201020000 
0001000100000002 
0102010000000000 
0000000000000000 
045E755142400101 
0101010000000000 
0000000000000000 
0000070507000000 
00000000083444C4 
44400D7901112D70 
050D750509700000 
0D70000000384540 
4020014E322E3140 


624 octets 


sm 


35E 
680 
9FD 
D57 
OEA 
484 
837 
BBA 
F25 
2A2 
5F9 
93F 
c92 
FEO 
333 
69C 
9F6 
D3C 
098 
3F0 
718 
A2B 
D3B 
048 
35B 
66B 
978 
c8B 
F9B 
2AB 
5BB 
8cB 
BDB 
F06 
270 
5B9 
8F3 
co3 
F13 
223 
539 
84E 
B62 
E76 
186 
4D2 
7E5 
AF5 
E18 
156 
486 
800 
B43 
E97 


036: 
037: 
038: 
039: 
03A: 
038: 
03C: 
030: 
03E: 
03F: 
040: 
041: 
042: 
043: 
044: 
045: 
046: 
047: 
048: 
049: 
O4A: 
04B: 
O4C: 
04D: 
O4E: 
04F: 
050: 


CONNLEX 


000: 
001: 
002: 
003: 
004: 
005: 
006: 
007: 
008: 
009: 
004: 
008: 
00C: 
000: 
00€: 
00F: 
010: 
011: 
012: 
013: 
014: 
015: 
016: 
017: 
018: 
019: 


084E794142400000 
00000000002E4559 
3200000000000000 
0000000000000026 
5556587008365556 
5810083645464830 
0832414248700024 
5655587008345655 
5810083446454830 
0C3042414C700024 
5556587008355654 
5810083546444830 
0C3142404C700025 
5455587008355455 
5810083544454830 
0C3140414C700875 
1414187000044972 
40000E3159454E30 
OC7A0F7949400024 
5554587000084A71 
40000C523A262D10 
0424587458400875 
1415187000094A70 
4000083544454830 
0c3140414C300074 
5655545000054071 
40000 


0123456789ABCDEF 


34F4E4E4C4548502 
802E000000000000 
4B200C5192900000 
F020000000000000 
04D000DD00AB000D 
934F4E4053519974 
25140535291FF17F 
8FDF8E0840948508 
508F2EA21AF68DCD 
CE017F8FBC631590 
AF2D6018D91FB08D 
07450306A0E 10A8F 
D9630781011A90A2 
18FDC6305F011A49 
D8D271308DC5E203 
1358D393904BFFF6 
BFFF85360103AFFF 
SAFFF8438F871F03 
2760AB58F155F040 
C16A146132EA1301 
63309A85AF215CF1 
6FAOD55F8FDF8E08 
42958508527FFE84 
1860508511087DEE 
1097C0F10A750F10 
B1128FBBCE011086 


1E7 
525 
83A 
B52 
EB1 

202 
543 
8A0 
BEF 
F44 
2A1 

5F0 
946 
CAO 
FEE 
350 
6A1 

A01 

D79 
0D5 
436 
78D 
ADD 
E21 
189 
4EO 
5D9 


344 octets 


sm 


392 
éc1 
A16 
D3E 
OC5 
443 
7c2 
B65 
F33 
2E4 
6A2 
A21 
D96 
130 
4CB 
876 
c47 
005 
38F 
706 
AA2 
E80 
214 
594 
913 
cB3 


O1A: 
018: 
01C: 
01D: 
O1E: 
01F: 
020: 
021: 
022: 
023: 
024: 
025: 
026: 
027: 
028: 
029: 
02A: 
028: 
O2C: 
02: 


STRUC2 


000: 
001: 
002: 
003: 
004: 
005: 
006: 
007: 
008: 
009: 
O0A: 
008: 
00C: 
000: 
O0E: 
00F: 
010: 
011: 
012: 
013: 
014: 
015: 
016: 
017: 
018: 
019: 
O1A: 
018: 
01C: 
01D: 
O1E: 
01F: 
020: 


11187AS5E0BF88414 
5S0A7A11BB8FBBCE01 
007DCEAF7AFABFBB 
CE0103112AFBE68F 
BBCE011186011B7A 
5EOBF8840450A7A1 
O1AFBAFACCAB881C 
81C3030E01B0D909 
40E4136E2134AC3B 
47550A&7A0D59F11 
111487011B7A5E0B 
F8850450A7A101DB 
8FBBCE0110870C08 
618186370871E0B7 
ASBOBF8450A7A11B 
8FB7CE086250B648 
2281C832D1A47580 
B47160A3F93F288D 
84A8015640E4F154 
448DF 


0123456789ABCDEF 


3545255534230202 
802E000000000000 
FOCO001E246448B00 
FB30087000000000 
072600D230F6600D 
61019600D520A960 
0D90069500D554E4 
44249C4541465546 
4B25540554144544 
955E44594C454975 
8494C454341FF303 
0014040202020202 
CC230BF354727573 
64757275602D4963 
7D6164713686CFF0 
38DE6A208DA2C207 
SFFFE1E34420FE1E 
06820FE1E16710FE 
1E36410006351185 
012165002218580C 
0208D82D20185011 
371FAA8F2AC21554 
135719FFE1E46DDF 
007E408FB394031A 
8966A18F8EC20AA6 
BB6BB679AF6B1072 
207430313F966D08 
F8EC2070208FEA23 
044B8D271308F9DF 
30870935CE8F9DF3 
0870A230187340E6 
8161371FAA8F2153 
494811135942008D 
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04B 
408 
803 
BB5 
F4E 
2E7 
écc 
A4C 
DDB 
167 
4F2 
887 
C2D 
FA3 
351 
6F9 
A74 
ETA 
198 
2c7 


1541 octets 


sm 


342 
671 
9EE 
D38 
OA6 
408 
784 
AEA 
E44 
1cB 
S4E 
871 
BFF 
F6C 
306 
6AE 
A7E 
E19 
16A 
4B6 
814 
BBC 
F85 
30A 
éc1 
A7E 
DEC 
196 
53A 
8E9 
C50 
FDC 
33B 








021: 53E201554135018F 6A1 057: 08F7C1B120AC0482 36F 08D: 361B198F21441340 44F 


022: D96307B9F7BCE4F9 A73 058: 860013028160E4E6 6D6 08E: 11371F1C6F214513 7C2 
023: 000068401818FB39 DDB 059: B20871B13018160E A45 08F: 5011F1C6F21471E9 B4C 
024: 40310F9628131049 134 05A: 4E681014A14F9E29 DD7 090: 95F1458D681F0D01 EEO 
025: 62F03112962606F1 49D 058: E3048160E4EACBOE 18B 091: 4A136C213401614E 24D 
026: 0626F1371F088F21 81C 05C: 4694C201B698F214 50E 092: 80CF208F5341080D 5D7 
027: 57413594E40018D5 B8D 05D: 61340163DFF47BFF 8BF 093: F80C0201F088F215 95F 
028: 30308D0745014B7B EFF 05E: 18776232873336D8 C2D 094: 501F765F214787D3 CE7 
029: D4474D61708F6242 284 O5F: 1EFCFF9FAFF70D28 029 095: 1321088F1C811137 04D 
02A: 0510220300339464 5BE 060: FE3F804033098169 3B0 096: C2D731FE8F99A805 40A 
02B: 23653037C4F4F405 93B 061: 4752DB13416314A8 728 097: CA795F16114A311E 7A1 
O2C: 2765203B3554C454 CA7 062: F3E3203471006010 A81 098: 962606DA016114A3 B14 
02D: 34452B6110397584 007 063: 449000066210C0C8 DDF 099: 124962023F340644 E6C 
02E: 94C45429749014B8 382 064: 0C020A87309A0B80 15B 09A: 16365426462F8F89 1EE 
02F: D303507D80390245 6EB 065: DO7A6276C26021D2 4E0 09B: 0B14F78F3E32034F 58A 
030: 8454E4296FDF7D54 A9E 066: CFFC9AFF14A7EFO4 8D4 O9C: 5006850441503664 8D3 
031: 5213754C43554276 E00 067: D0300902725D1313 C30 09D: 016B30001F088F21 C3F 
032: 8cCF153331A896631 183 068: 479E17F915E07032 FB7 O9E: 5708F62420850811 F9E 
033: 8F0E160754017252 4EE 069: 7CB16FDO08D84A803 368 09F: 760E01D01D60F707 324 
034: 27240313F966C137 85D 06A: 10677C175A168C05 6E5 0A0: FO5DO0A8058008028 698 
035: 0245F40227722017 BBO 06B: DBFFD3AFF7212776 ABB OA1: 6D002767000C0C80 A09 
036: 114B7F108FCE2503 F46 06C: 15E078F174816BCF E61 OA2: CF2070211F088F21 D8E 
037: 1C221770017166AF 2B9 06D: 207A4264A051CFF0 1F1 0A3: 554798E1817FCE63 131 
038: 8D324508D229506F 63D O6E: 1AFF78616FAF89BF 5D1 OA4: OF7A7E16514A7D1D 4DD 
039: EFF83EFF18770351 9FA 06F: FO0AFF1361B698F2 980 0A5: 51160D031547FD04 850 
03A: 678F871F030998A1 D8A 070: 144314473611B698 CEO 0A6: 716BFC314670D048 BDE 
03B: 331E0962828F83DB 120 071: F21461347C014E07 056 OA7: 02261906CA031267 F3F 
03C: 07C151B178F2137C 4A9 072: D9179216C406C6F7 3F0 0A8: BB043F215C70C0C8 2E0 
030: 2144164140AC25C0 812 073: CBFFAEAFF7/BEO4DE 80F 0A9: OCF20AC714A33462 66E 
03E: 7DF4AC2B461B098F BD3 074: 257AD164305ABFF6 BBA OAA: 6238F890B15C3312 9EC 
03F: 2154477C47015267 F34 075: 1BFF2376C1602020 F39 0AB: 4966911611520301 D33 
040: D1578C4165312414 29F 076: 33301E8DA93908D5 2CB DAC: 90242302902C1ACB OAE 
041: A962F016114A7D43 627 077: 045078411877E811 62E 0AD: 8126F0F31267F504 42D 
042: 5606D9031A8966A3 9AB 078: 36134D787DD13210 9A8 OAE: 702150275DD18178 7A5 
043: 1F088F2161156415 D18 079: 88F1C8111378B70C D3B OAF: 1E6F4E3146704048 B25 
044: 541607AA41361B08 08A O7A: C28B39B86F7786D2 OEE 0B0: EOCOC80CF2071BD1 ED2 
045: 8F2156413480DF77 412 07B: 71F0B7F22030214B 46B 0B1: 6514A745C4028129 240 
046: 605444B5778414A3 783 O7C: OEO6AOC4951FD55F 817 082: 0681789D18172504 5A7 
047: 13F962E0227B4058 B01 07D: 214713517E143F43 B80 083: 0077DD6B0E632C14 93D 
048: 24F3267F30551161 E6B 07E: 4412E08A68313610 EES 0B4: A96200011F495F21 CB2 
049: 7F54237F205C0432 1E6 07F: A1341FC65F2143E6 27E 085: 43E41418F156818B 035 
04A: 1617C4414A311F96 560 080: E6BBAC18F95EF011 641 OB6: E911F088F2157013 3B4 
04B: 6C016114A6D6F673 8ED 081: A1348FB7EF011213 9D4 087: 11C01550018DD449 726 
04C: F7BF365B280CF201 C98 082: 08D7E4708F871F08 D73 OB8: 01FF85F214717414 AAS 
04D: 361B698F21441B09 018 083: FD8DBO04ACO05D097 12C 0B9: 38A200cC14113115 EOA 
O4E: 8F21524A4C413109 38E 084: C00110AC097C0001 499 OBA: 701F088F21550011 16D 
O4F: 8F45DB0044B11298 71A 085: 3098161851328F31 7FA 0BB: E064600000FB3000 4CD 
050: ODF1298F534D0042 AA4 086: F808161341650110 B4D OBC: 00000000D1091BFF 83A 
051: 065B08D51DBO0AC78 E4C 087: 88FE3F8043430981 ED9 OBD: D8204C7FFD00030A BDA 
052: F83DB0137C21351B 1DF 088: 694783DB134AF015 269 OBE: FFD61035BFF9B005 F97 
053: 678F2146841840AC 564 089: A51657630118AE5A 5F1 OBF: 6BFFD73414355426 32B 
054: OB8AA218A8868A601 8F0 08A: F2AE9BF2BF2BF2BF 9EE OCO: 7540435544639464 697 
055: 8505418A8C25C38B C82 088: 233FE1E9720061AE D86 OC1: 367C4F4F40506B35 A27 
056: 670851D618414213 FDF O8C: 1B198F2146134011 O0E8 OC2: 54C4543445161FF D73 


JPC 59 Page 32 





Le Journal JPC est le bulletin de liaison entre les 
membres de l'Association "PPC Paris", régie par la loi 
de 1901. Le Club est éditeur de JPC, et son siège 
social est au 56, rue Jean-Jacques Rousseau, 75001 
Paris. 


La maquette de ce numéro a été préparée et réalisée 
par Pierre David, Janick Taillandier et Jacques 
Baudier grâce à un système comprenant un HP71B, 
un lecteur de disquettes HP9114A, un HP9807A, 
deux HP9154 et une imprimante LaserJet. 


Les dessins sont de Jean-Jacques Dhénin et Paul 
Courbis. 


Directeur de la publication : Pierre David 
Numéro ISSN : 0762 - 381X 


Veuillez adresser toute correspondance à : 
PPC Paris, BP 604, 75028 Paris Cedex 01. 


Imprimé par Copy-Express, 42 86 91 94. 


ENGLISH SUMMARY 
JPC 59 - NOVEMBER 1988 


We had the opportunity to see the new HP-42S to be available in a near future. It looks like an HP-27S but 
it has an enhanced HP-41 instruction set and 8 KBytes ram. Unfortunately, the 1/0 capabilities are limited 
to infra-red printer. 


On page 3, Eric Gengoux describes changes to his railway ticket price calculation program. See JPC 36 and 
JPC 44 for more details. He also presents an HP-27 expression for the SOLVE menu. 


The HP-28 columns begins with correction to the article about curve fitting published in JPC 57. Next, 
Dang-Trung Khoi-Nguyen gives us a method to build custom characters and as an example, a well known 
logo. Last, Serge Vaudenay gives us on page 7 his version of the HP-28C (1CC) ASSEMBLEUR program, as well 
as à function à similar to the Forth word. 


HP-41 columns includes a matrix and system solving programs. Two versions are provided, one using 
synthetic programming, the other one standard programming, but slightly slower. Thanks to Jean-Marc 
Baiïllard. 


For HP-75 and VisiCalc users, Eric Gengoux provides 5 extension functions in order to process 
chronological series. 


The Hp-71 pages begin with a new version of the Connett graph plotter program from Guy Toublanc. This 
version is 16 times faster than the previous one from Stefano Tendon (from Italy), thanks to a Lex file 
equivalent to the old internal loop. 


And last, Bruno Gïil plots a function on the HP-71 LCD display and allows you to use it as a window on the 
full graph. It uses only MathRom and the STRUC2 Lex. 


Until next month, 


Happy Programming and JPC reading ! 


